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Abstract 

The  energetic  ionic  liquids  formed  by  the  1,2,4-triazolium  cation  family  and  dinitramide 
anion  are  investigated  by  ab  initio  quantum  chemistry  calculations,  to  address  the 
following  questions:  How  does  substitution  at  the  triazolium  ring’s  nitrogen  atoms  affect 
its  heat  of  formation,  and  its  charge  delocalization?  What  kind  of  ion  dimer  structures 
might  exist?  And,  do  deprotonation  reactions  occur,  as  a  possible  first  step  in  the 
decomposition  of  these  materials? 
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Introduction 


Recent  years  have  seen  a  great  enthusiasm  for  ionic  liquids  (ILs)  in  regard  to  their  utility 
as  novel  and  “green”  solvents  '  ,  since  their  ability  to  serve  as  new  solvent  environments 
is  quite  remarkable.  IL  properties  include  high  electrical  conductivities,  high  dielectric 
constants,  non-existent  vapor  pressures,  and  often  excellent  thermal  and  chemical 
stabilities3  and  therefore  offer  the  prospect  of  hosting  new  catalytic  behavior4'5.  A  short 
list  of  IL’s  novel  solvent  properties  would  include  different  IL  systems  that  are  either 
completely  miscible  or  immiscible  with  water4,  the  ability  to  form  emulsions6,  and  tuning 
of  the  dielectric  constant  by  absorption  of  large  mole  fractions  of  supercritical  CO?7. 

Two  recent  computer  simulations  of  l-butyl,3-methylimidazolium  (bmim  )/PFis"  show  at 
least  two  different  time  scales  for  diffusion  in  such  solvents8,  and  yield  information  on 
preferred  structures9.  An  experimental10  study  of  the  interaction  of  water  with  the  surface 
layer  of  this  IL,  known  to  be  immiscible  with  water,  observed  structural  changes 
consistent  with  the  latter  simulation. 

In  spite  of  this  easily  understandable  focus  on  the  novel  solvent  properties  of  stable  ILs,  it 
is  good  to  remember1 112  that  ILs  are  neither  new  materials,  nor  required  to  be  stable,  nor 
limited  to  use  as  solvents.  Indeed,  the  well-known  and  undesired  “red  oil”  appearing 
during  Friedel-Crafts  reactions  is  an  ionic  liquid  containing  an  AI2CI7"  anion.  A  12°C 
melting  point  for  CH3CH2NH3  NO3"  was  reported  as  early  as  1914  ,  and  a  dark  red- 
brown  oil  Et3NH+CuCL’  was  discovered  in  1963 14.  King13  and  Wilkes16  and  co-workers 
investigated  systems  such  as  LiCl/AlCC  (m.p.  145°C)  starting  in  the  1960s,  with  the  goal 
of  developing  molten  liquid,  high  current  batteries.  Electrochemical  properties  of  the 
1, 4-dimethyl- 1, 2, 4-triazolium  chloride/AlCC  system  have  been  reported17.  A  key 
breakthrough  came  with  Wilkes’  design  using  MNDO  theory  of  the  1-ethyl, 3- 
methylimidazolium  (emim  )  cation  in  1982  ,  achieving  a  25°C  melting  point  for 
emimVAlCLf.  A  decade  later  the  Wilkes  group19  investigated  changing  the  anion  to  NO3" 
,  BF4',  and  acetate,  with  the  latter  freezing  to  a  glass  at  -45°C.  This  work  helped 
established  the  current  IL  design  principle  that  the  very  large  lattice  energy  of  a 
conventional  salt  such  as  NaCl  (m.p.  801°C)  can  be  lowered  to  near  or  even  below  room 
temperature  by  using  cations  and  anions  differing  in  size,  by  delocalizing  the  charge 
within  each  ion  to  reduce  point  charge  forces,  and  by  exploiting  asymmetrically  shaped 
ions  to  entropically  disrupt  the  lattice3.  Many  ILs  used  as  solvents  are  derived  from 
Wilkes’  emim1  in  that  they  contain  a  substituted  imidazolium  ring,  for  example  the 
increasingly  popular  bmim  .  The  recent  development  of  the  bis(triflyl)amide  anion20  has 
made  it  easy  to  create  numerous  IL  systems  that  avoid  crystallization  altogether,  to  form 
glasses  at  temperatures  well  below  0°C.  Very  recently”  it  has  been  shown  that  amino 
acids  can  serve  as  the  counterion  to  emim+,  demonstrating  that  the  anion  need  not  be 
inorganic.  Two  additional  non-solvent  applications  of  ILs  can  be  mentioned.  One  is  the 
development  of  an  ionic  liquid  crystal  system22.  Electrospray  of  the  stable  IL  emim+/BF4' 
has  been  proposed  for  use  as  a  highly  controllable  microthruster  '  . 
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Room  temperature  ionic  liquids  are  normally  considered  chemically  and  thermally  stable, 
since  they  are  often  based  on  the  imidazolium  or  other  stable  cations,  and  may  also 
contain  nearly  inert  perfluoro  anions.  However,  it  is  of  interest  to  consider  chemically 
unstable  ionic  liquids26,  which  may  have  interesting  applications.  In  particular,  this  work 
will  consider  the  use  of  both  more  energetic  cations  and  anions.  The  cations  considered 
are  the  triazoliums,  five  membered  rings  containing  three  nitrogens  rather  than  the  two 
found  in  imidazolium.  Note  that  a  dinitrosubstituted  triazole  was  recently  used  as  a 
triazolate  anion  in  an  energetic  low  melting  salt  but  herein  we  consider  only  triazolium 
cations.  The  anion  considered  replaces  the  stable  ions  used  in  solvents  such  as  BF4'  or 
PF6-  with  the  energetic  dinitramide  ion  ,  NfNCFV. 

Accordingly,  the  first  topic  addressed  in  the  current  paper  is  an  attempt  to  assess  charge 
delocalization  in  triazolium  ions.  Second,  the  heats  of  formation  of  various  substituted 
triazolium  cations  are  predicted.  Next,  by  considering  ionic  pairs  (one  cation  and  one 
anion),  possible  structure  motifs  in  the  bulk  ionic  liquids  are  probed.  Finally,  as  a  first 
attempt  to  consider  the  possible  decomposition  reactions  of  an  energetic  IL,  we  report  on 
the  deprotonation  reactions.  Relevant  literature  citations  for  each  topic  will  be  given 
below,  as  each  is  considered  in  turn.  Related  research  on  tetrazolium  based  energetic 
ionic  liquids  will  be  presented  elsewhere  . 

Method  of  Calculation 

As  all  species  considered  here  are  closed  shell  in  nature,  the  majority  of  the  calculations 
are  done  at  the  closed  shell  Hartree-Fock  (RHF)  or  second  order  perturbation  theory 
(MP2)  levels.  To  facilitate  n  bonding  analysis,  some  multi-configuration  (MCSCF) 
calculations  were  carried  out  to  access  the  inherent  interpretability  of  this  wavefunction 
in  terms  of  localized  active  orbitals32,  although  the  RHF  determinant  in  all  cases  remains 
at  least  92%  of  the  MCSCF  function.  Mulliken  charges33  and  charges  fitted  to 
electrostatic  potentials34  were  not  found  to  yield  useful  ring  charge  delocalization 
information.  The  basis  set  employed  for  all  calculations  is  6-3 1++G(d,p),  using  diffuse 
functions  on  all  atoms  (including  the  cationic  species)  so  that  these  functions  are  present 
for  the  anions.  Some  G2  thermochemistry  calculations36  are  performed  using  the 
standard  methods  and  basis  sets  employed  in  the  construction  of  that  model.  All 
computations  were  perfonned  with  the  GAMESS  electronic  structure  package  ’  ,  except 
the  G2  calculations,  which  employed  the  Gaussian94  code  . 

Charge  Delocalization  in  Triazolium  Cations 

Triazolium  based  ionic  liquids  are  typically  prepared  by  quatemization  of  a  triazole 
(often,  by  protonation  of  a  triazole  ring  nitrogen  by  a  strong  acid),  followed  by  a 
metathesis  reaction  to  replace  the  anion  of  the  generated  triazolium  salt  by  anions  more 
suitable  to  fonning  ionic  liquids.  The  triazole  or  triazolium  can  be  1,2,3  or  1,2,4  (see 
below  for  the  IUPAC  ring  numbering),  but  only  the  latter  is  considered  in  the  present 
paper.  The  synthesis  of  4-amino- 1,2, 4-triazole  (4AT)  has  been  known  since  194440. 
Reactions  which  transform  the  NH2  functionality  at  N4  into  various  substitutions  at  Ni  are 
well  known,  see  for  example41.  The  IR  spectra  of  the  parent  1,3-imidazole  H4C3N2, 
1,2,4-triazole  H3C2N3,  and  even  tetrazole  H2CN4  were  published  in  195642,  developing  an 
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IR  fingerprint  for  protonation  at  ring  nitrogens  (in  1,3  imidazolium  cation).  Raman 
spectra  of  chlorhydrate  salts  of  triazolium  bromide43  give  similar  fingerprints  for 
triazolium  NH  ring  vibrations.  Recent  experimental  work44  on  the  quaternization  of  4AT 
shows  that  protonation  occurs  at  Ni,  and  conversely,  quaternization  of  1 -amino- 1,2,4- 
triazole  leads  to  protonation  at  N4.  An  energetic  salt  (m.p.  95°C)  of  4AT  with  nitroform 
ion  (CXNCfiV)  was  reported  in  196645,  and  triazolium  salts  with  the  energetic  perchlorate 
CIO4"  were  given  in  196946.  The  group  of  Pevzner47  has  worked  on  the  synthesis  of 
energetic  triazoles  and  related  compounds  for  two  decades.  More  recently  Shreeve  and 
coworkers48'53  have  reported  the  synthesis  of  numerous  ionic  liquids  containing 
imidazolium,  triazolium,  and  tetrazolium  cations.  A  number  of  crystal  structures  of  the 
bromide  and  nitrate  salts  of  substituted  triazoliums  are  now  available54. 

Substitution  in  triazolium  is  possible  at  both  nitrogen  and  carbon  atoms,  with  the  known 
substituents  at  the  N  atoms  being  primarily  H,  alkyl,  or  NH2  (in  any  combination  at  N 1 
and  N4),  although  the  interesting  nitroimide  (-N'-N02)  group  has  been  placed  at  N455. 
Routes  to  disubstitution  at  Ni  and  N2  rather  than  the  more  typical  Ni  and  N4  are  also 
known56.  Energetic  substituents  at  carbon  include  both  nitro  (NO2)47’57  and  azido  (N3) 
groups.  The  azido  group  was  introduced  into  an  azole  in  190  547  and  such  azoles  were 
recently  quatemized  into  ionic  liquids  .  It  is  also  possible  to  introduce  halogen  at  one  or 
both  carbons  in  the  neutral  triazoles,  see  for  example  .  The  following  “back  of  the 
envelope”  estimation  predicts  that  substitution  at  nitrogen  should  give  a  higher  heat  of 
formation  than  carbon  substitition.  A  textbook  tabulation59  of  typical  bond  strengths 
shows  that  breaking  a  CH  bond  should  cost  about  25  kJ/mol  more  than  breaking  an  NH 
bond.  However,  the  payback  for  forming  a  C-R  bond  compared  to  an  N-R  bond  more 
than  recovers  this  cost.  C-R  bond  energies  exceed  N-R  bond  energies  by  202,  157,  58,  or 
41  kJ/mol,  respectively,  for  R=C,  O,  N,  F  single  bonds.  Thus  carbon  substituted 
triazoliums  would  be  expected  to  have  lower  heats  of  formation,  and  therefore  only 
nitrogen  substitution  is  considered  here.  Quantitative  tests  of  this  assumption  in  the 
tetrazolium  system  will  be  reported  elsewhere29. 

Obviously,  triazolium  based  energetic  ionic  liquids  are  experimentally 
realistic26’47’52’53’54’60,  so  it  is  relevant  to  ask  how  substitution  on  the  triazolium  may  affect 
their  properties.  The  two  questions  about  the  triazolium  cation  that  are  addressed  first  are 
the  charge  delocalization  in  nitrogen  substituted  1,2,4-triazolium  cations,  and  in  the  next 
section,  their  gas  phase  thermodynamic  stabilities.  Recall  that  greater  charge 
delocalization  is  considered  to  be  a  key  factor  in  producing  lower  melting  points.  The 
triazolium  species  considered  here  include  the  more  interesting  1,4-substituted 
compounds,  but  also  include  the  1,2-substituted  forms  for  completeness, 
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Typically  only  one  substituent  R,  R’,  or  R”  is  considered  here,  with  the  other  ligand 
usually  being  H.  The  substituents  R  are  chosen  to  include  the  simple  series  H,  F,  OH, 
NH2,  and  CH3,  and  more  energetic  groups  such  as  CN,  N3,  NF2,  N02,  CH=CH2,  or 
OCH3.  Table  1  gives  the  relative  energies  for  monosubstitution.  Note  that  substitution  at 
the  1  position  of  1,4  isomers  is  always  energetically  favored  over  substitution  at  the  4 
position,  although  sometimes  by  only  a  narrow  margin.  NMR  experiments44  show  that 
for  R=NH2,  the  only  observed  cation  is  the  1,4  isomer,  R’=H.  Because  the  1,2  isomers 
are  always  higher  in  energy,  they  will  not  be  considered  further  (except  for  the  inclusion 
of  a  few  R=R”=H  results). 


Table  1.  MP2/6-31++G(d,p) 

electronic  energy  differences, 

in 

kj/mol . 

H 

F 

OH  NH2 

CH3  CN 

n3 

nf2 

no2 

c2h3 

och3 

1,4 

isomer  (R'=H) 

R=  0 

0 

0  0 

0  0 

0 

0 

0 

0 

0 

4,1 

isomer  (R=H) 

R'=  0 

16 

20  21 

12  0.2 

22 

7 

16 

12 

22 

1,2 

isomer  (R"=H) 

R=  5  6 

73 

65  67 

55  56 

63 

58 

57 

53 

64 

The  geometry  of  the  triazolium  ring  is  remarkably  constant.  Compared  to  the  parent  1 ,4- 
substituted  triazolium,  the  other  1,4  and  4,1  substituted  systems  change  the  five  ring  bond 
lengths  by:  Ni-N2=  -0.029  to  0.014,  N2-C3=  -0.005  to  0.009,  C3-N4=  -0.008  to  0.016,  N4- 
Cs=  -0.008  to  0.016,  and  Cs-Ni=  -0.004  to  0.01 1  A,  at  the  MP2/6-31++G(d,p)  level. 

It  is  now  convenient  to  introduce  the  resonance  forms  that  will  used  to  discuss  the  charge 
distribution  in  1,4-substituted  triazolium  ions,  and  in  the  next  section,  for  constructing 
isodesmic  reactions61.  In  order  to  assess  the  relative  importance  of  these  resonance 
structures,  a  charge  and  bonding  analysis  of  the  parents  (R=R’=  H): 
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w 
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will  be  presented.  Resonance  structures  I  and  II  differ  in  the  position  of  double  bonds, 
and  in  the  location  of  the  cationic  charge.  The  important  design  criterion  for  ionic  liquids 
that  they  should  exhibit  substantial  charge  delocalization  is  therefore  the  question  as  to 
the  relative  contribution  these  two  resonance  forms  make  to  the  molecular  wavefunction. 
The  1 ,2-substituted  triazolium  has  two  equivalent  resonance  structures  III  whereas  the 
neutral  triazoles  IV  and  V  produced  by  deprotonation  have  only  the  single  uncharged 
Lewis  diagrams  shown.  All  four  species  are  6n  e‘  systems. 
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Two  routine  procedures  for  assigning  atomic  charges  were  tested.  Chemical  intuition 
strongly  suggests  that  any  resonance  between  I  and  II  implies  substantial  positive  charge 
on  both  Ni  and  N4.  However,  Mulliken  population  analysis33  of  either  RHF  or  MP2 
densities  showed  large  changes  at  the  ring  atom  as  a  function  of  substituent,  for  example 
the  charge  at  Ni  for  R=CN  is  -0.33  while  for  R=N3  it  is  +0.27.  Since  no  chemically 
meaningful  patterns  can  be  seen,  the  data  are  omitted.  Secondly,  fitting  the  MP2 
electrostatic  potential  with  atomic  point  charges34  also  did  not  give  sensible  results.  For 
example,  for  ten  of  the  eleven  cases  for  R  substitution,  the  charges  at  Ni  are  0.30  to  0.55, 
while  N4  actually  takes  on  a  negative  value  ranging  -0.21  to  -0.50.  For  R’  substitution, 
the  charges  reverse,  with  Ni  quite  close  to  zero  in  ten  of  1 1  compounds,  while  N4  ranges 
0.12  to  0.27.  Since  both  of  these  nitrogens  should  be  positive,  these  results  are  regarded 
as  reflecting  the  procedure  of  fitting  to  the  total  electrostatic  potential,  with  each  atom’s 
charge  merely  a  fitting  parameter. 

However,  the  relative  importance  of  the  resonance  forms  can  be  assessed  by  analysis  of 
the  MCSCF  wavefunctions,  in  which  the  six  n  electrons  are  considered  to  occupy  five 
active  n  orbitals.  The  closed  shell  determinant  comprises  92-93%  of  all  such  MCSCF 
wavefunctions  considered,  but  the  MCSCF  wavefunctions  afford  a  much  more  revealing 
analysis  of  what  the  n  electrons  are  doing  than  ordinary  Hartree-Fock.  The  analysis 
procedure32  begins  by  localizing  the  five  active  MCSCF  orbitals,  resulting  in  five  atom 
centered  p  orbitals,  then  regenerates  the  MCSCF  density  matrix  for  these  atomic-like 
MCSCF  MOs.  Diagonal  elements  of  the  localized  orbital  density  matrix  may  be 
considered  to  be  the  number  of  electrons  occupying  each  particular  p  orbital.  Off- 
diagonal  elements  are  n  bond  orders,  with  positive  and  negative  values  reflecting  bonding 
or  antibonding,  and  with  the  magnitude  of  the  number  giving  the  strength  of  the 
interaction.  Table  2  contains  the  results  of  such  an  analysis  for  the  1,4  and  1,2 
dihydrogen  triazolium  ions,  and  their  deprotonated  forms.  The  energetics  of 
deprotonation  is  addressed  later,  as  the  present  focus  is  an  assessment  of  the  resonance 
stabilization. 
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Table  2.  MCSCF/ 6-31++G (d, p)  71  orbital  analysis. 

71  orbital  populations  total  MP2 


Ni 

n2 

c3 

n4 

c5 

bonding 

Energy 

1 , 4  ion 

I/II 

1 .43 

1.15 

0.98 

1.49 

0.95 

1  .  91 

-241 . 932210 

1 , 2  ion 

III 

1.50 

1.50 

0.94 

1 . 12 

0.94 

1  .  94 

-241 . 910862 

neutral 

IV 

1 . 14 

1 . 14 

1.05 

1 . 62 

1.05 

1 .  98 

-241 . 574345 

neutral 

V 

1.59 

1.20 

1.03 

1.16 

1 . 02 

1 .  95 

-241.585265 

adjacent 

(1,2) 

K  bonding 

next  neighbor 

(1, 

3)  71  antibonding 

Ni-N2 

n2-c3 

c3-n4 

N4-C5 

Cs-N, 

Ni-C3 

n2-n4 

C 

3-C5  n4-n4 

c5-n2 

1 , 4  ion 

I/II 

0 . 42 

0.78 

0.48 

0.58 

0  .  67 

-0 .14 

-0 .19 

-0 

.18  -0.28 

-0.23 

1 , 2  ion 

III 

0.39 

0  .  62 

0  .  63 

0  .  63 

0  .  62 

-0.30 

-0.20 

+  0 

.05  -0.20 

-0.30 

neutral 

IV 

0 . 44 

0 .77 

0.46 

0.46 

0 .77 

-0.10 

-0 .24 

-0 

.24  -0.24 

-0.10 

neutral 

V 

0 . 42 

0.73 

0.54 

0.70 

0.53 

-0.21 

-0 . 15 

-0 

.05  -0.25 

-0.31 

If  resonance  structure  I  dominated  over  II,  Ni  would  possess  one  n  electron,  and  N4  two  n 
electrons.  Conversely,  structure  II  assigns  two  n  electrons  to  Ni  and  one  to  N4  Table  2 
shows  that  N 1  and  N4  possess  1 .43  and  1 .49  n  electrons,  respectively,  implying  nearly 
equal  weights  for  resonance  forms  I  and  II.  This  means  the  fonnal  plus  charge  in  I  and  II 
is  roughly  equally  shared  by  these  two  nitrogens.  Similarly,  Ni  and  N2  have  1.50  n 
electrons  in  the  1-2  isomer,  where  symmetry  implies  the  resonance  forms  III  must  have 
exactly  equal  contribution. 

The  bonding  analysis  in  the  second  part  of  Table  2  supports  the  implication  of  nearly 
equal  resonance  weights  obtained  from  the  n  populations.  None  of  the  adjacent  n 
bonding  interactions  is  as  small  as  0  or  as  large  as  1,  as  implied  in  the  drawings  I-V, 
when  any  particular  n  bond  is  absent  or  present.  Focusing  on  N4C5  and  C5N1,  which 
exchange  a  n  bond  between  resonance  forms  I  and  II,  we  see  that  neutral  molecule  IV  has 
a  smaller  n  bond  order  (0.46)  between  N4C5,  while  neutral  V  has  a  more  complete  bond 
here  (0.70),  with  1 ,4-triazolium  being  intermediate  at  0.58.  Similarly,  the  0.67  bond 
order  for  C5N1  in  1,4  triazolium  is  also  intermediate  between  the  0.77  and  0.53  for  the 
neutrals. 

Every  next  nearest  neighbor  interaction  is  slightly  antibonding  (negative  interaction), 
with  but  one  exception.  This  is  a  common  characteristic  in  aromatic  systems,  for 
example,  all  meta  interactions  in  various  benzene  analogs  are  found  to  be  antibonding62, 
and  this  merely  reflects  the  inherent  delocalization.  If  we  simply  sum  all  five  bonding 
interactions  plus  all  five  antibonding  interactions,  the  resulting  “total  bonding” 
interaction  reported  in  Table  2  is  remarkably  close  to  the  2  n  bonds  drawn  in  each 
diagram  I-V. 

Table  3  reports  the  n  populations  and  total  bonding  of  various  1,4  substituted  triazolium 
cations.  It  may  be  seen  that  apart  from  the  ligands  CN  and  NF2,  there  are  virtually 
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constant  populations  at  Ni  and  N4  although  the  values  are  somewhat  changed  from  those 
of  the  parent  (R=R’=H),  Ni=1.43  and  N4=T.49,  reported  in  Table  2. 


Table  3.  MCSCF/ 6-31++G (d, p)  n  orbital  analysis  of  substituted  triazolium 
cations . 

n  orbital  populations 


R 

R' 

Ni 

N 

2 

c3 

N 

u 

C 

5 

total  bonding 

F 

H 

1  . 

37 

1 . 

16 

99 

1  . 

49 

1 . 

00 

1 . 

91 

OH 

H 

1  . 

32 

1 . 

17 

99 

1  . 

52 

1 . 

00 

1 . 

92 

nh2 

H 

1 

.36 

1 

.  15 

1 

.00 

1  . 

.51 

.  97 

1 

.  93 

ch3 

H 

1 

.33 

1 

.15 

1 

.00 

1  . 

.53 

1 

.00 

1 

.  92 

CN 

H 

1  . 

54 

1 . 

13 

1  . 

00 

1  . 

40 

94 

1 . 

90 

n3 

H 

1 

.32 

1 

.  17 

.  99 

1  . 

.52 

1 

.00 

1 

.  90 

nf2 

H 

1 

.48 

1 

.13 

1 

.00 

1  . 

.44 

.  95 

1 

.  92 

no2 

H 

1 

.33 

1 

.15 

1 

.01 

1  . 

.52 

1 

.00 

1 

.  91 

C2H3 

H 

1 

.33 

1 

.15 

1 

.00 

1  . 

.52 

.  99 

1 

.  92 

och3 

H 

1 

.31 

1 

.16 

1 

.00 

1  . 

.53 

1 

.00 

1 

.  93 

H  F 

1  . 

44 

1 . 

14 

1  . 

01 

1  . 

43 

99 

1 . 

91 

H  OH 

1  . 

44 

1 . 

14 

1  . 

01 

1  . 

42 

98 

1 . 

92 

H  NH2 

1 

.44 

1 

.  14 

1 

.02 

1 

.43 

.  97 

1 

.  93 

H  CH3 

1 

.46 

1 

.  14 

1 

.02 

1 

.39 

.  99 

1 

.  93 

H  CN 

1  . 

34 

1 . 

13 

98 

1  . 

60 

95 

1 . 

92 

H  N3 

1 

.43 

1 

.  14 

1 

.01 

1 

.45 

.  97 

1 

.  92 

H  NF2 

1 

.37 

1 

.13 

.  99 

1 

.56 

.  95 

1 

.  93 

H  N02 

1 

.46 

1 

.15 

1 

.01 

1 

.40 

.  99 

1 

.  90 

H  C2H3 

1 

.45 

1 

.  14 

1 

.01 

1 

.41 

.  99 

1 

.  92 

H  OCH3 

1 

.46 

1 

.  14 

1 

.02 

1 

.39 

.  99 

1 

.  92 

CN  NH2 

1 

.57 

1 

.13 

1 

.01 

1 

.35 

.  94 

1 

.  91 

nh2 

CN 

1 

.30 

1 

.  14 

0 

.  98 

1 

.  62 

.  95 

1 

.  91 

Taking  CN  rather  than  NF2  from  the  two  substituents  that  affect  n  populations,  and  taking 
NFF  as  typical  of  the  others,  we  include  two  doubly  substituted  systems  as  the  final 
entries  in  Table  3.  Indeed,  these  two  molecules  seem  to  represent  the  limit  to  which  we 
can  push  n  electron  localization  onto  Ni  or  onto  N4.  Of  course  low  melting  points  are 
associated  with  charge  delocalization,  so  the  only  prediction  that  can  be  made  regarding 
lower  melting  points  is  to  use  ligands  other  than  NF2  and  CN.  The  slightly  more  stable  R 
substituted  molecules  appear  to  be  a  bit  less  charge  delocalized  than  the  R’  substituted 
systems.  The  total  n  bond  orders  of  all  molecules  remain  almost  unchanged,  at  nearly 
two  full  n  bonds  in  the  rings. 

Thermochemistry  of  Triazolium  Cations 

Speaking  of  energetic  cations  implies,  of  course,  that  their  heats  of  formation  are  high. 
The  gas  phase  heats  of  formation  will  be  presented  for  both  Ni  and  N4  monosubstituted 
triazolium  ions,  as  this  quantity  is  independent  of  the  nature  of  the  counter  anion.  The 
most  efficient  computational  procedure  for  finding  accurate  heats  of  formation  is  based 
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on  isodesmic61  (bond  type  conserving)  reactions.  In  the  present  case,  since  we  have  two 
resonance  structures  I  and  II  given  above,  we  must  consider  two  isodesmic  reactions: 


I  +  2  CH4  +  5  NH3  ->  2  H3C-NH2  +  2  H2C=NH  +  H2N-NH2  +  H2N-R’  +  H3N-R+  AHi 

II  +  2  CH4  +  5  NH3  ->  2  H3C-NH2  +  2  H2C=NH  +  H2N-NH2  +  H2N-R  +  H3N-R’+  AHn 

which  differ  by  exchanging  R  and  R’  in  the  amine  and  protonated  amine  products.  By 
having  the  same  number  of  CN  single  and  double  bonds,  and  NN,  CH,  and  NH  single 
bonds  on  each  side  of  these  isodesmic  reactions,  one  hopes  for  improved  accuracy  (due  to 
error  cancellations)  in  the  computed  enthalpy  of  reaction,  obtained  in  the  present  case  by 
MP2/6-3 1++G(d,p).  Combining  the  total  reaction  energy  with  the  experimental  heats  of 
formation  of  all  smaller  molecules  in  these  two  reactions  allows  extraction  of  a 
theoretical  heat  of  formation  for  both  I  and  II.  It  will  be  shown  that  these  two  heats  are 
always  very  similar,  and  their  average  can  therefore  be  taken  as  a  heat  of  formation  of  the 
actual  resonance  hybrid  of  I/II. 

Since  many  of  the  R  groups  considered  here  are  rather  unstable,  many  of  the  amine  and 
protonated  amine  heats  of  formation  are  not  available  from  experiment.  These  are 
therefore  computed  by  the  G2  procedure,  and  compared  to  available  experimental 
information,  in  Table  4. 


Table  4.  G2  energies,  proton  affinities,  and  heats  of  formation  of 
amines  and  protonated  amines,  at  298°K.  G2  total  energies  are  in 
Hartree,  while  all  other  energies  are  in  kj/mol.  PA  is  the  proton 
affinity  at  the  amine  nitrogen,  whereas  PA'  is  the  proton  affinity  at 
the  preferred  protonation  site  (when  different) .  The  G2  energies  for  CH4 
and  H2C=NH  are  -40.407060  and  -94.459413,  respectively. 

==  RNH2  ==  ==  RNH3+  == 


R 

G2 (RNH2) 

G2  (RNH3+) 

PA 

PA' 

PA' 

AHf 

AHf 

AH£ 

AHf 

G2 

exp . 

G2 

G2 

exp . 

G2 

exp 

H 

-56 . 454841 

-56.777598 

847 . 4 

853.6 

-47 . 6 

-45.9 

626.4 

630 

F 

-155.529266 

-155 . 817447 

756.6 

-33.2 

732 . 1 

OH 

-131 . 529103 

-131 . 837695 

810.2 

-47 . 3 

-42 

664 . 4 

nh2 

-ill  .  676216 

-112 . 003327 

858 . 8 

853.2 

96.7 

95 . 4 

759.7 

770 

ch3 

-95 . 662525 

-96.003321 

894 . 8 

899.0 

-25 . 4 

-23 . 0 

601.7 

611 

CN 

-148 . 554808 

-148 . 817548 

689.8 

805.6 

810 . 8 

135.7 

134 

967 . 8 

n3 

-219.783477 

-220 . 087659 

798.6 

418.8 

1142 . 1 

nf2 

-309.883456 

-310.162220 

731 . 9 

-27 . 8 

762.2 

no2 

-260 .702675 

-260 . 976511 

719.0 

757 . 4 

746.6 

-5 . 8 

797.1 

c2h3 

-133 . 686373 

-134 . 012498 

856.2 

898 . 9 

917 . 7 

55 . 9 

29 

721.5 

och3 

-170.740391 

-171.060178 

839.6 

844.8 

-34 . 6 

-26 

647 . 7 

The  G2  energies  including  thennal  corrections  to  298°K  are  obtained  in  the  standard 
way36.  Experimental  proton  affinities63  agree  with  G2  values  to  within  the  2  kcal/mol 
accuracy  expected  for  the  G2  method  (that  is,  to  about  10  kJ/mol),  except  for  the  case 
R=C2H3.  Note  that  three  of  the  R  groups  contain  more  favorable  protonation  sites  than 
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the  amine  nitrogen.  Heats  of  formation  for  all  amines  and  protonated  amines  are 
obtained  from  the  G2  energies  using  G2  atomic  energies363,  and  G2  heats  of  formation  at 
0°K  of  atoms366,  and  the  standard  G2  thermal  corrections  for  atoms366.  A  value  of 
1528.085  kJ/mol  was  used  for  the  heat  of  formation  of  the  gas  phase  proton64,  with  a 
thennal  correction  of  25.93  kJ/mol  (5RT/2).  Experimental  heats  of  formation65  shown  in 
Table  4  also  agree  well  with  the  results  from  G2  theory,  again  with  the  exception  of  the 
vinyl  substituent  C2H3. 

The  computed  heats  of  reaction  for  the  two  isodesmic  reactions  above  are  found  by 
combining  the  MP2/6-31++G(d,p)  energy  difference  for  the  reaction,  the  scaled  zero 
point  energies  and  other  thermal  factors  (ideal  gas,  rigid  rotor,  harmonic  oscillator 
partition  functions),  to  yield  a  theoretical  AH°(298)  for  reactions  1  and  2.  Given  the  G2 
derived  heats  of  formation  presented  in  Table  4,  it  is  easy  to  extract  the  heats  of 
formation  of  the  resonance  forms  I  and  II,  and  then  their  average  value,  which  are 
presented  in  Table  5. 


Table  5.  Heats  of  formation  in  kJ/mol  at  298  for  gas  phase 
monosubstituted  1 , 2 , 4-triazolium  cations,  at  298°K.  I  and  II  refer  to 
particular  resonance  hybrids,  with  their  average  values  being  the 
predicted  heats.  The  average  is  given  in  both  kJ/mol  and  kJ/g. 


For  R"=H  cations  For  R=H  cations 


R 

AHf°(I) 

AHf°  (II) 

ave 

rage 

R" 

AHf°(I) 

AH£°  (II) 

average 

H 

822.2 

822.2 

822 

11.7 

H 

822.2 

822.2 

822  11.7 

F 

895.2 

904 . 8 

900 

10.2 

F 

920 . 4 

910 . 8 

916  10.4 

OH 

809.4 

822.1 

816 

9.5 

OH 

842 . 9 

830.2 

837 

9.7 

nh2 

898 . 0 

905 . 4 

902 

10 . 6 

nh2 

927 . 0 

919.6 

923  10.9 

ch3 

770 . 6 

773 . 7 

772 

9.2 

ch3 

786.2 

783.1 

785  9.3 

CN 

1108 . 3 

1109 . 1 

1109 

11.7 

CN 

1109.6 

1108.8 

1109  11.7 

n3 

1272 . 8 

1288 . 8 

1281 

11.5 

n3 

1310 . 9 

1294.9 

1303  11.7 

nf2 

919.5 

940.4 

930 

7 . 7 

nf2 

947 . 4 

926.5 

937 

7 . 7 

no2 

971.9 

935 . 6 

954 

8.3 

no2 

951.7 

987 . 9 

970 

8 . 4 

c2h3 

878 . 9 

884 . 4 

882 

9.2 

c2h3 

898 . 1 

892 . 6 

895  9.3 

och3 

798.4 

809.9 

804 

8 . 0 

och3 

831 . 7 

820.2 

826 

8.2 

The  only  comparison  with  experiment  that  is  possible  is  for  the  parent,  R=R”=H,  for 
which  the  experimental  heat  of  formation  is  835  kJ/mol65.  The  closeness  of  the  present 
computed  result,  822  kJ/mol,  to  this  value  is  encouraging.  Note  that  very  similar  heats  of 
formation  are  derived  from  the  isodesmic  reaction  of  each  resonance  structure,  with  NO2 
having  the  largest  variation  at  36  kJ/mol.  This  suggests  that  the  averaged  heats  of 
formation  given  in  Table  5  may  be  taken  as  an  estimate  of  the  heat  of  formation,  in  the 
gas  phase,  of  the  monosubstituted  triazolium  cations.  As  already  noted  with  regard  to  the 
0°K  energy  differences  given  in  Table  1,  substitution  at  N4  is  slightly  higher  in  energy 
than  at  Np  While  all  of  the  energies  in  Table  5  are  large  positive  numbers,  because  they 
are  gas  phase  reference  states,  on  a  per  weight  basis  the  parent  and  the  two  very  energetic 
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groups  N3,  CN  stand  out.  Note,  however,  that  the  experimentally  accessible  protonated 
4AT  (R=NH2)  is  close  behind  this  trio.  Perhaps  the  most  realistic  synthetic  target  for 
high  energy  content  suggested  by  Table  5  is  CN  substitution  at  nitrogen,  since  N3 
substitution  would  require  creating  a  ring  exo  NN  single  bond.  The  rest  of  this  paper 
concentrates  on  two  cations,  namely  the  parent  and  NH2  substitution  at  N4,  because  of 
their  high  energy  content  and  their  experimental  existence. 


Ion  Deprotonation  and  Protonation  Energetics 

In  the  next  section  we  will  consider  pairs  of  ions,  either  the  parent  1  -H,4-H- 1,2,4- 
triazolium  or  l-H,4-NH2-l,2,4-triazolium  cation  combined  with  a  dinitramide  anion,  and 
the  proton  transfer  reactions  between  these.  Here,  the  energy  costs  for  deprotonating 
either  ring  hydrogen  from  an  isolated  triazolium,  and  the  energy  gains  for  protonating  an 
isolated  dinitramide  anion  at  all  possible  sites  are  presented,  as  these  factors  combine  to 
determine  the  energetics  of  most  of  the  ion  and  neutral  pair  geometries  located. 

The  energies  required  to  deprotonate  1-H,4-H  or  1-H,2-H  substituted  triazolium  cations, 
to  produce  all  possible  neutral  products  can  be  found  in  Table  6.  These  energies  were 
obtained  by  combining  the  MP2  total  energies  of  the  ions  and  neutrals  given  in  Table  2 
with  MP2  zero  point  energies  and  thermal  corrections.  The  tabulated  experimental 
value63  of  the  proton  affinity  of  1-H- 1,2, 4-triazole  is  in  reasonably  good  agreement  with 
the  computed  value,  but  experiment  does  not  give  the  preferred  protonation  site.  The 
primary  literature  cited  by  the  standard  table  used  MP2/6-3  lG(d,p)  energies  at 
RHF/STO-3G  geometries  to  interpret  their  experimental  result,  finding  an  energy 
preference  of  58  kJ/mol  for  protonation  at  N4  over  N2,  nearly  identical  to  the  55  kJ/mol 
reported  in  Table  6.  A  much  older  paper43  using  CNDO/2  electrostatic  potentials  leads  to 
the  same  conclusion,  that  1-H- 1,2 ,4-triazole  preferentially  protonates  at  N4  rather  than 
N2.  However  a  recent  paper67  using  B3LYP/6-3 1 1+G(3df,3pd)  at  B3LYP/6-3  lG(d,p) 
structures  reported  a  theoretical  gas-phase  basicity  attributed  to  protonation  at  N2,  which 
is  quite  possibly  a  typographical  error68. 
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Table  6.  Cation  deprotonation  and  anion  protonation  energetics 
(kj/mol).  Values  are  AH°(298°K),  in  kj/mol. 


deprotonate  at  Nil 
deprotonate  at  N4: 
deprotonate  at  N2 : 


1 , 4-triazolium  -> 
1 , 4-triazolium  -> 
1 , 2-triazolium  -> 


4-triazole  +  H 
1-triazole  +  H 
1-triazole  +  H 


927 

901  (exp  .  =88 6hl) 
846 


protonate  at  central  N:  N(N02)2  +  H+  ->  HN(N02)2  -1312 

protonate  at  terminal  0:  N(N02)2~  +  H+  ->  N(N0(0H))N02  -1293 


Deprotonation  from  Ni  in  1-H,4-NH2- 1,2, 4-triazolium  cation  requires  941  kj/mol, 
compared  to  the  927  kj/mol  required  for  the  parent,  so  N4  amine  substitution  might  affect 
the  hydrogen  bonding  ability  of  the  more  tightly  bound  Ni  hydrogen. 

Table  6  also  contains  data  for  the  protonation  of  the  energetic  ion  NINCkV,  named 
dinitramide28.  This  ion  is  considered  here  because  in  addition  to  its  high  energy  content, 
it  has  several  sites  for  the  negative  charge  to  delocalize,  and  therefore  triazolium  and  4- 
amino-triazolium  cations  form  a  low  melting  salt  and  an  oil,  respectively,  when  combined 
with  dinitramide  anion"  .  It  is  known  from  theoretical  calculations"  that  the  free 
dinitramide  ion  has  a  C2  geometry,  but  Raman  spectra  show  that  it  is  found  in 
dinitramide  salts  without  any  local  symmetry.  Protonation  of  dinitramide  has  previously 
been  considered69,  and  our  present  results  are  essentially  the  same.  Protonation  at  the 
central  nitrogen,  giving  a  structure  without  any  symmetry,  is  favored  by  19  kj/mol  over 
protonation  at  the  terminal  oxygens.  There  are  four  isomers  for  protonation  on  the 
tenninal  oxygens,  all  of  which  have  Cs  symmetry,  with  the  NO2H  group  coplanar  with 
the  plane  of  the  three  nitrogens.  The  most  stable  of  these  has  the  H  atom  on  the  oxygen 
atom  trans  to  the  other  NO2  group,  which  is  perpendicular  to  the  N3  plane, 

M»0 

On^n'"0 


1 

H 

so  the  energy  of  this  isomer  was  used  in  Table  6.  Trial  structures  with  the  proton  placed 
near  the  N  of  an  NO2  group  optimized  towards  protonation  at  the  oxygens.  Both  the 
dinitramide  ion  and  its  protonated  forms  are  found  to  be  rather  flexible,  being  able  to 
rotate  NO2  groups  and  deform  various  bend  angles. 

Note  that  the  protonation  energy  of  dinitramide  is  larger  than  the  deprotonation  energy  of 
triazolium  cations,  indicating  that  the  anion  may  be  able  to  pull  a  proton  away  from  the 
cation,  a  theme  that  will  be  developed  in  the  next  section. 
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Structures  of  Triazolium  and  Dinitramide  ion  pairs 

Now  consider  the  possible  dimer  geometries  for  a  single  triazolium  cation  and  a  single 
dinitramide  anion.  This  is  closely  related  to  simulations  of  both  small  clusters  and  bulk 
materials  by  the  Thompson  group  for  the  explosive  materials  ammonium  nitrate70'72, 
ammonium  dinitramide  '  ,  and  hydroxylammonium  nitrate  .  This  research  ,  and 
earlier  computations  on  ammonium  dinitramide  in  the  Morokuma  group  ,  show  that  one 
NH4+  cation  and  one  dinitramide  anion  spontaneously  transfer  a  proton  to  produce  neutral 
products.  Although  a  minimum  for  an  ionic  dimer  pair  exists  at  the  RHF  level  ,  both 
MP278  and  DFT75  computations  find  that  only  a  hydrogen  bonded,  neutral  pair 
(NH3)(HN(N02)2)  exists.  However,  this  is  due  to  the  small  size  of  the  ion  dimer  cluster, 
as  the  ion  tetramer  (NH4+)2(N(N02)2’)2  or  even  the  ammonia  adduct  to  the  ion  dimer 
(NH4+)(N(N02)2')NH3  are  minima  on  the  correlated  potential  energy  surfaces,  due  to 
increased  numbers  of  strong  Coulombic  interactions75. 

Although  there  is  no  question  that  bulk  materials  will  consist  of  ionic  species,  these  dimer 
proton  transfer  reactions  and  the  neutral  pairs  that  result  are  likely  to  be  a  primary 
decomposition  mechanism,  since  decomposition  probably  occurs  simultaneously  with 
vaporization.  Alavi  and  Thompson  have  studied  the  decomposition  of  protonated 
dinitramide,  HN(N02)2,  and  related  systems  such  as  its  parent  anion  and  ammonium 
complexes,  and  give  references  to  earlier,  relevant  experimental  and  theoretical  work. 
They  find  that  the  two  experimentally  observed  channels  for  ammonium  dinitramide 
decomposition,  N20  or  N02  products,  each  have  an  activation  energy  of  about  35 
kcal/mol,  but  that  N20  product  can  be  formed  from  the  solid  state  ionic  structures, 
whereas  the  NO?  channel  opens  only  after  vaporization  leads  to  proton  transfer  to  neutral 
precursors. 

In  addition  to  their  relevance  to  the  decomposition  of  energetic  ionic  liquids,  dimer 
structures  may  provide  some  insight  into  the  bulk  liquid  structure.  Here  we  consider  in 
detail  the  structures  of  l-H,4-H-l,2,4-triazolium  cation  with  dinitramide  anion,  and  less 
thoroughly,  l-H,4-NH2-l,2,4-triazolium  cation  (quaternized  4AT)  with  dinitramide. 
Experimentally26,  the  former  system  is  found  to  be  a  low  melting  salt  (m.p.  75°C),  while 
the  latter  is  a  yellow  oil  (m.p.  20°C).  Unfortunately  no  crystal  structure  of  the  former 
system  is  available  for  comparison  to  the  present  results,  but  X-ray  structures  of  similar 
materials  suggest  two  possible  types  of  ion  interactions:  coplanar  or  three  dimensional. 
Hydrogen  bonding  is  the  most  important  dimer  interaction,  and  this  will  be  maximized 
when  the  H-bond  is  linear.  Examples  of  crystal  structures  showing  coplanar  interactions 
include  the  salts  1 -H,3-N3,4-H-l, 2, 4-triazolium  nitrate52  or  l-H,4-H-l,2,4-triazolium 
chloride'  .  Of  course,  solids  contain  more  than  just  simple  dimer  interactions,  and  many 
body  electrostatic  forces  may  lead  to  more  three-dimensional  structures,  in  which  as 
much  H-bonding  as  possible  is  retained,  while  allowing  ions  to  be  surrounded  by  more 
oppositely  charged  ions.  An  example  of  this  is  l-H,4-H-l,2,4-triazolium  perchlorate26 
which  contains  some  nearly  linear  NH-0  interactions,  but  many  more  nearly 
perpendicular  NH-0  and  CH-0  interactions,  resulting  in  a  3D  lattice.  A  set  of  crystal 
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structures54  for  doubly  substituted  triazoliums  (organic  groups  at  Ni  and  amine  at  N4)  as 
bromide  or  nitrate  salts  have  predominantly  3D  organizations,  controlled  by  hydrogen 
bonding  between  the  amine  and  the  anion. 

As  anticipated  in  the  previous  section,  the  lowest  energy  neutral  (proton  transferred) 
dimer  is  a  hydrogen  bonded  system  combining  the  lowest  energy  triazole  with  the  lowest 
energy  protonated  dinitramide,  shown  in  the  lower  right  corner  of  Figure  1.  The 
hydrogen  bond/dipole  interaction  is  fairly  strong,  so  separating  these  neutral  species  from 
each  other  requires  5 1  kJ/mol.  It  is  possible  to  locate  ionic  dimer  structures  at  the 
Hartree-Fock  level,  and  the  corresponding  ion  dimer  in  Figure  1  lies  19  kJ/mol  above  the 
proton  transferred  neutral  dimer.  Of  course,  separating  the  ions  to  infinity  requires  a 
much  larger  energy  than  separating  neutrals,  422  kJ/mol  relative  to  the  neutral  pair.  The 
similar  magnitude  of  electrostatic  stabilization  that  would  result  from  bringing  additional 
ions  together  (such  as  Alavi  and  Thompson75  found  for  the  ion  tetramer 
(NFLt+)2(N(N02)2  )i)  means  that  for  triazolium  based  clusters  with  more  than  two  ions, 
ionic  structures  can  be  presumed  to  lie  below  neutral  structures.  Figure  1  also  contains  an 
estimate  of  the  charge  separation,  found  by  summing  the  Mulliken  charges  over  all  atoms 
in  each  ion  at  the  RHF  level,  showing  that  the  ionic  structure  is  indeed  quite  ionic.  MP2 
single  point  calculations  at  the  RHF  geometry  have  a  somewhat  smaller  Mulliken  charge 
on  the  ions,  ±0.85  e'. 

Figures  2,  4,  and  5  show  that  the  triazolium/dinitramide  system  is  very  rich  in  both  ionic 
and  neutral  structures.  Figure  2  contains  all  ionic  geometries  found,  which  can  transfer 
protons  through  the  transition  states  shown  in  Figure  3,  to  yield  the  neutral  structures 
given  in  Figure  4.  The  position  of  each  dimer  in  Figures  2,  3,  and  4  is  kept  consistent,  so 
that  any  individual  proton  transfer  reaction  may  be  easily  followed.  The  common  zero  of 
energy  for  all  structures  in  Figures  2-5  is  the  most  stable  neutral  dimer  (already 
considered  in  detail  in  Figure  1),  which  can  be  found  at  the  top  center  of  Figure  4.  Note 
that  the  orientation  of  the  ring  is  kept  the  same  in  all  geometries  shown,  so  that  the 
reader’s  eye  can  focus  on  the  changing  position  of  the  dinitramide  partner  around  the 
ring. 

The  top  half  of  Figures  2  and  4  involves  hydrogen  bonding  at  the  N4  position  of  the  ring, 
where  the  deprotonation  of  triazolium  is  easiest.  The  bottom  half  of  these  Figures  gives 
H-bonded  structures  at  Nj,  where  deprotonation  is  more  difficult.  Thus  all  energies  in  the 
bottom  half  of  Figure  4  are  higher  than  the  top  half.  Within  either  half,  the  energy 
ordering  can  be  rationalized  by  noting  that  the  central  N  of  dinitramide  is  the  protonation 
site  preferred  over  the  tenninal  oxygens. 

Figure  5  collects  a  few  additional  neutral  structures  in  which  the  H  atom  in  the  protonated 
dinitramide  is  not  hydrogen  bonded  to  the  triazole.  Three  of  the  five  structures  are  non- 
planar,  and  in  fact,  the  stacked  structure  at  the  bottom  right  is  found  to  exist  only  at  the 
MP2  level.  The  bottom  left  figure  has  the  central  nitrogen  of  the  diimide  nearly  coplanar 
with  the  ring,  so  its  hydrogen  is  not  participating  in  a  hydrogen  bond.  Most  of  these  have 
a  rather  high  energy,  and  will  not  be  considered  further. 
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The  minima  presented  in  Figures  2  and  4  are  mainly  “coplanar”,  meaning  that  the  three 
nitrogen  atoms  in  the  dinitramide  moiety  lie  near  the  plane  of  the  rings,  although  they  are 
never  precisely  coplanar  (structures  have  Ci,  not  Cs  symmetry).  These  structures  were 
located  by  geometry  optimization  from  40  trial  geometries,  generated  by  placing  the  5 
possible  pairs  (dinitramide  ion  with  triazolium  ion,  or  either  O  or  N  protonated 
dinitramide  with  either  possible  triazole)  oriented  with  the  three  nitrogens  perpendicular 
to  the  ring,  in  all  8  octants,  about  7A  from  center  to  center.  These  40  starting  geometries 
were  thus  biased  against  the  predominantly  coplanar  structures  that  result  from  geometry 
optimization,  which  are  therefore  clearly  preferred  by  the  hydrogen  bonding  contacts. 
After  consideration  of  the  results  of  these  “random”  geometry  searches,  additional,  more 
rational  searches  were  carried  out.  A  total  of  65  geometry  searches  for  ionic  or  neutral 
pairs  were  performed,  resulting  in  the  29  minima  presented  in  Figures  2,  4,  and  5.  It  is 
unlikely  that  all  minima  have  been  found!  However,  only  one  three  dimensional  ion 
cluster  and  matching  neutral  complex  was  found,  shown  at  the  bottom  right  corner  of 
Figures  2  and  4.  One  additional  non-planar  minimum  is  not  shown,  namely  a  minimum 
similar  to  the  27  kJ/mol  structure  in  Figure  4,  in  which  the  protonated  dinitramide  is 
rotated  to  be  perpendicular  to  the  triazole  ring,  at  a  relative  energy  of  3 1  (30  kJ/mol  at  the 
MP2  level).  This  structure  was  omitted  since  it  appears  to  be  just  an  internal  rotational 
form  of  the  planar  structure  depicted,  as  no  such  perpendicular  ionic  structure  could  be 
optimized. 

MP2  energy  computations  at  the  RHF  geometries  give  the  second  set  of  energy  values 
given  in  Figures  2-5.  Although  4  of  the  1 1  ionic  geometries  given  in  Figure  2  have  lower 
energies  at  the  RHF  level  than  their  neutralized  partners  in  Figure  4,  MP2  computations 
preferentially  stabilize  the  neutral  forms,  so  all  1 1  ionic  geometries  lie  higher  in  energy 
than  their  neutral  fonns.  It  may  be  seen  by  comparing  Figure  2  to  Figure  3  that  although 
there  are  small  barriers  at  the  Hartree-Fock  level  for  the  proton  transfer  reactions  that 
convert  ionic  pairs  to  neutral  pairs  (for  example,  a  barrier  of  8  kJ/mol  for  the  most  stable 
neutral  dimer,  featured  in  Figure  1),  in  all  cases,  the  MP2  energies  of  the  ionic  fonns  in 
Figure  2  are  slightly  higher  than  that  of  the  transition  states  in  Figure  3.  Therefore,  nearly 
all  structures  given  in  Figure  2  are  just  shoulders  on  the  MP2  potential  surface,  and  MP2 
geometry  optimization  leads  to  spontaneous  proton  transfers  to  structures  like  those 
shown  in  Figure  4.  However,  the  two  ionic  molecules  in  the  bottom  row  of  Figure  2  do 
possess  shallow  minima  at  the  MP2  level.  The  third  set  of  relative  energies,  where  given, 
are  MP2  energies  at  MP2  structures.  Note  that  the  MP2  barriers  for  proton  transfers  from 
the  two  stable  ionic  pairs  are  very  small,  just  0.6  and  2.8  kJ/mol. 

Figure  6  contains  the  only  dimer  structure  found  for  the  1,2  substituted  triazolium  ion. 
Figure  6  is  actually  part  of  the  same  potential  energy  surface  as  the  structures  in  Figures 
2-5,  and  is  singled  out  only  because  it  contains  a  1,2  substituted  triazolium.  Both  neutral 
and  ionic  structures  are  clearly  energetically  competitive,  lying  below  many  of  the 
structures  given  in  Figures  2  through  5.  The  5 1  kJ/mol  (see  Figure  1)  needed  to  separate 
the  two  monomers  of  the  most  stable  neutral  pair,  followed  by  reattachment  at  N2,  is  an 
upper  limit  on  the  energy  required  to  access  the  neutral  structure  in  Figure  6.  Although 
the  1 ,2  substituted  triazolium  is  more  easily  deprotonated  than  1 ,4  substituted  triazolium 
(see  Table  6),  an  MP2  ionic  minimum  does  exist,  because  the  dinitramide  is  able  to  H- 
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bond  to  both  protons  rather  than  just  one  in  the  1,4  ionic  complexes.  Clearly  other 
neutral  structures  with  hydrogen  bonds  between  the  ring  N2  and  oxygen  atoms  in  the 
dinitramide  moiety  may  exist,  at  somewhat  higher  energy,  but  these  were  not  searched 
for. 

The  available  crystal  structures54  for  bromide  or  nitrate  salts  of  disubstituted  triazoliums 
(Ni  organic,  N4=NH2)  show  three  dimensional  structures  resulting  from  H-bonding  by 
the  NH2  groups.  Accordingly,  the  results  of  a  far  less  extensive  search  for  ionic 
structures  that  might  exist  in  the  yellow  oil,  l-H,4-NH2-l,2,4-triazolium  dinitramide  is 
given  here.  The  top  part  of  Figure  7  shows  the  result  of  eight  “random”  searches  as 
outlined  above,  with  the  two  ions  started  about  7 A  apart,  in  all  8  octants.  These  searches 
did  not  lead  to  any  minima  involving  hydrogen  bonding  at  N 1  although  clearly  this  search 
procedure  is  not  exhaustive.  Instead,  the  minima  that  were  found  show  that  hydrogen 
bonding  at  the  amine  group  is  quite  favorable,  accompanied  by  additional  interactions 
between  the  ring  and  dinitramide.  The  existence  of  six  such  similar  minima  illustrates 
the  large  degree  of  flexibility  in  dinitramide  ion.  All  six  structures  survive  MP2 
geometry  optimization,  which  is  not  suprising  since  deprotonation  at  the  amine  would 
produce  a  zwitterion,  not  a  neutral  triazole.  An  explicit  search  for  a  structure  with  a 
dinitramide  ion  hydrogen  bonded  at  its  central  nitrogen  to  the  Ni  ring  proton  resulted  in 
the  location  of  a  ionic  minimum  at  the  RHF  level,  shown  in  the  lower  left  corner  of 
Figure  7,  at  an  energy  comparable  to  the  NFF  bound  structures.  The  previous  section 
reported  a  slightly  higher  deprotonation  energy  for  this  ring  proton  (941  vs.  927  kJ/mol  in 
the  parent),  and  the  RHF  barrier  seems  more  substantial  than  in  the  parent.  Nonetheless, 
the  MP2  energy  at  the  RHF  saddle  point  is  already  1 1  kJ/mol  below  the  MP2  single  point 
energy  of  the  ionic  structure,  and  thus  geometry  optimization  of  this  ionic  structure  leads 
to  spontaneous  proton  transfer  to  the  neutral  pair  in  the  lower  right  of  Figure  7,  whose 
energy  is  substantially  below  the  NH2  bound  ionic  pairs.  In  short,  the  chemistry  of  the 
ring  proton  is  qualitatively  the  same  as  that  of  the  parent  triazolium,  while  the 
introduction  of  the  amine  changes  the  chemistry  at  N4. 

When  combined  with  dinitramide  anion,  the  unsubstituted  l-H,4-H-l,2,4-triazolium 
cation  has  a  preference  for  coplanar  geometries,  whereas  l-H,4-NH2-l,2,4-triazolium 
adopts  more  three  dimensional  shapes.  The  amine  group  allows  less  ordered  structures  to 
exist,  giving  a  plausible  explanation  for  why  the  former  system  forms  a  low  melting  salt 
(m.p.  75°C)  but  the  latter  is  a  yellow  oil  (m.p.  20°C).  Of  course,  an  X-ray  structure  for 
the  1-H,4-H-1,2,4  triazolium  dinitramide  salt  would  be  welcome,  to  confirm  whether  the 
dimer  geometries  reported  here  are  relevant  in  the  bulk  solid. 

Summary 

The  results  presented  in  this  paper  address  several  different  aspects  of  energetic  ionic 
liquids.  The  conclusions  reached  include  the  following: 

1 .  The  contribution  of  the  two  resonance  structures  in  the  unsubstituted  triazolium  cation 
is  roughly  equal,  with  approximately  3/2  n  electrons  found  on  both  charged  nitrogens. 
Substitution  at  ring  nitrogen  atoms  produces  little  change  in  the  extent  of  this  charge 
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delocalization,  as  shown  by  both  very  small  geometry  changes  in  the  rings,  and  by 
MCSCF  analysis  of  the  resonance  hybrids.  The  two  substituents  that  cause  the  most 
degradation  to  resonance  are  CN  and  NF2,  but  to  only  a  small  extent.  It  appears  that 
substitution  will  not  hinder  the  ability  of  these  cations  to  form  ionic  liquids. 

2.  A  number  of  energetic  and  less  energetic  substituents  at  the  ring  nitrogens  were 
investigated,  and  in  all  cases  1 ,4  substituted  triazoliums  are  lower  in  energy  than  1 ,2 
substituted  triazoliums.  There  is  a  small  energy  preference  for  monosubstitution  at  Ni 
over  N4  in  the  1,4  substituted  triazoliums.  The  heats  of  fonnation  for  substituted  gas 
phase  triazolium  cations  were  presented,  and  the  CN  substituent  is  suggested  as  a 
synthetic  target,  over  N3,  since  the  former  does  not  require  creating  a  ring  exo  NN  bond, 
but  rather  a  NC  bond. 

3.  Triazolium  dinitramide  systems  are  rich  in  structure.  The  numerous  geometries  and 
their  relative  energies  are  understandable  by  consideration  of  hydrogen  bonds  and 
protonation  energetics.  The  presence  of  only  small  energy  barriers,  or  often,  spontaneous 
transfer  of  protons  from  ionic  dimers  to  produce  neutral  pairs  implies  that  deprotonation 
is  a  fundamental  mechanism  for  triazolium  decomposition. 

The  most  important  question  that  was  not  addressed  in  the  present  paper  relates  to  the 
decomposition  pathways  of  these  energetic  ionic  liquids,  beyond  suggesting 
deprotonation  as  an  initiating  step.  Little  seems  to  be  known  about  the  decomposition  of 
triazoles  or  triazolium  ions,  so  experiments  that  detect  products  from  their 
decomposition,  or  further  theoretical  work  would  be  welcome.  For  example,  would  a 
triazole  ring  produced  by  deprotonation  subsequently  extrude  molecular  nitrogen?  In 
contrast,  a  considerable  amount  of  information  about  the  decomposition  of  the 
dinitramide  ion  is  available79. 
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Figure  1.  Relative  energies  (beneath  each  structure)  of  the 
most  stable  neutral  dimer,  and  its  corresponding  ionic 
dimer.  Structures,  energies,  and  Mulliken  charges  are 
from  RHF/6-3 1++G(d,p).  As  discussed  in  the  text,  the 
ionic  dimer  does  not  survive  geometry  optimization  at  the 
MP2  level. 


Figure  2.  Ionic  pairs  in  the  parent  1,2,4-triazolium  and 
dinitramide  system,  located  at  the  RHF/6-31++G(d,p)  level. 
The  first  and  second  energy  values  are  from  RHF/6- 
31++G(d,p)  and  MP2/6-31++G(d,p)  at  RHF  geometries,  in 
kJ/mol.  The  third  value,  if  given,  is  at  the  fully  optimized 
MP2/6-31++G(d,p)  geometry.  The  zero  of  energy  for 
Figures  2-6  is  a  neutral  complex  presented  in  Figure  4. 
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Figure  3.  Transition  states  in  the  1,2,4-triazolium  and  dinitramide 
system,  connecting  ionic  pairs  to  neutral  complexes.  The  three 
energy  values  have  the  same  meaning  as  defined  in  Figure  2. 
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Figure  4.  Neutral  dimers  in  the  1,2,4-triazolium  and 
dinitramide  system,  resulting  from  proton  transfers.  The  three 
energy  values  have  the  same  meaning  as  defined  in  Figure  2. 
The  lowest  energy  structure  found  is  at  the  top  center,  and  is 
taken  as  the  energy  zero  for  Figures  2-6. 
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Figure  5.  Additional  neutral  complexes,  in  which  the 
dinitramide's  H  is  not  hydrogen  bonded  to  the  ring.  The  three 
energy  values  have  the  same  meaning  as  defined  in  Figure  2. 


(( 


Figure  6:  Structures  for  1-2  substituted  triazolium 
and  dinitramide:  ionic  mininum,  transition 
state, and  neutral  minimum.  The  three  energy 
values  have  the  same  meaning  as  defined  in 
Figure  2,  and  the  same  common  energy  zero  as 
defined  in  Figure  4. 
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Figure  7.  Structures  for  1 -hydro, 4-amino- 1,2, 4  triazolium  dinitramide. 
All  but  two  of  the  structures  shown  are  geometry  optimized  at  the 
MP2/6-31G++G(d,p)  level,  and  possess  positive  definite  hessians  at 
the  RHF/6-31++G(d,p)  level.  The  upper  and  lower  energies  are  from 
RHF  and  MP2,  respectively,  in  kJ/mol.  The  ionic  structure  and 
transition  state  in  the  bottom  row  are  RHF/6-31++G(d)  geometries. 
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